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ABSTRACT 

Three methods for generating outcomes on multivariate 
normal random vectors are presented. A comparison is made 
to determine which method requires the least computer 
execution time and memory space when utilizing the 
IBM 360/67. All methods use as a basis a standard Gaussian 
random number generator. Results of the comparison study 
indicate that the method based on triangular factorization 
of the covariance matrix generally requires less memory 
space and computer time than the other two methods. 
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I. INTRODUCTION 



Frequently it is desired to generate a sample of vectors 
derived from a specified multivariate normal distribution. 
For example, this need arises in simulation studies in which 
the model contains stochastic variates that are distributed 
according to a multivariate normal. Since a computer is 
frequently used to produce the sample vectors, the memory 
space and computer time requirements become important cost 
factors when considering various methods for "generating" 
random vectors. 

The problem of "generating" random vectors was discussed 
in 19^8 by Herman Wold [Ref. 1] . Wold describes a method 
for construction of samples from a multivariate normal 
distribution. Wold used a method based on triangularlzatlon 
of the covariance matrix. In 1962, Ernest M. Scheuer and 
David S. Stoller wrote a paper on the generation of multi- 
variate normal random vectors [Ref. 2] . They considered a 
method based on matrix equations and a method based on con- 
ditional distributions. The first method uses a technique 
similar to that of Wold's. Their conclusion was that the 
mathematics involved in the first method was simpler than 
that of the second method so the generation of random 
vectors is best accomplished by the matrix equation 
technique . 

The objective of this paper is to examine the above tv;o 
methods, hereafter called "triangularlzatlon" and 
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"conditioning" respectively, and to discuss a third method 
which is referred to as the "rotation" method. Computer 
programs for each method are presented and a comparison of 
memory space requirements and execution times on the 
IBM 360/67 computer is made for the three methods. The 
format for making this comparison study is similar to that 
used by M. E. Muller in his paper on univariate normal 
generators [Ref. 31 • 

The following notation is used in this study: X denotes 

a random variable whereas x is a real variable. A vector 
la denoted by a bar under the letter. Listed below are two 
examples of the above notation: 




In the discussion of the three methods, generators for 
N(0_,E) distributions are considered, where Z is a general 
(positive definite symmetric) matrix. There is no loss of 
generality in assuming the mean vector to be the zero 
vector since a random vector Y distributed N(u,E) may be 
generated by first generating X distributed N(^,E) and then 
taking Y = X + u. 
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II . METHODS 



A. ROTATION METHOD 

Given a positive definite symmetric p x p matrix Z, the 
function 

f(x;Z) = (I) 

is a p-variate normal density. A random vector X having 
this distribution has mean 0_ and variance covariance matrix 
Z [Ref. . Denote the distribution law corresponding to 
the density function (1) by N(OjZ). If X (with p components) 
is distributed N(£jZ), then Y = CX is distributed N(0_jCZC')j 
if C is nonsingular [Ref. 4]. The development of the 
rotation method is based on this fact. If a p x p matrix C 
is found such that CZC* yields a diagonal matrix, then 
elements on the diagonal of CZC’ are the variances of the 
corresponding components of Y. Since the covariance values 
(off diagonal elements of CZC’) are zero, the components of 
y are independent. 

For a symmetric matrix Z, there exists an orthogonal 
matrix C such that CZC’ = D where D contains the eigen 
values of Z on its diagonal [Ref. 51 • Associated with the 
eigen values are eigen vectors which make up the rows of 
the matrix C. Since Y = CX is N(£,CZC), and there exists 
a C such that CZC’ is the diagonal matrix D, it follows 
that the components of Y are independent normal variates 
with zero mean and variance values as given by the 
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corresponding eigen values of I. These variates may be 
generated using a univariate normal random number generator. 
A method of generating the random vector Y Is thus easily 
established. To obtain a generator for the random vector 
X, a transformation Is made using X = C~^Y. Since C is 
an orthogonal matrix, this simplifies to the equation 
X = C’Y. Standard matrix multiplication of the orthogonal 
matrix C’ and sample vectors generated on the random 
vector Y yields the desired random samples of the vector X. 

B. CONDITIONING METHOD 

A multivariate probability density function f may be 
written as; 

f(x^,...,x^) = f(x^|x 2 ,. . . ,x^) f(x 2 |x 2 ,. . . ,x^) 

••• ( 2 ) 

where f(x.|x....) denotes a conditional density function. 

If X^ can be generated from f(x^), then may be gener- 
ated by conditioning on X^. then be generated by 

conditioning on X^ and In a similar manner, all the 

components of X may be generated. In what follows, this 
general approach is applied to multivariate normal 
distributions . 

Given that the vector X is partitioned into two sub- 
vectors X^^^ and X^^^ and that the covariance matrix of X 
is correspondingly partioned into ^12’ ^22’ 

it can be shown [Ref h] that the conditional distribution 
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„ ^( 1 ) . „( 2 ) ( 2 ) . 1 .^u ^ ^ -1 ( 2 ) 

of X given X = x is normal with mean ^22 — 

and covariance matrix - E^^ ^22^ ^21* Several defini- 

tions are Introduced at this time as they are required in 
what follows. 

Let 



k,k 



k,k+l 



k ,n 



k+l,k k+l,k+l 



(k) _ 



k+1 ,n 



n ,k 



n,k+l 



n,n 



and let e[^^ denote the cofactor of o. . in E^^^: similarly 
E. . denotes the cofactor of a. , in E. Consider now the mean 
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and variance of the distribution determined by 
f (Xj^ I . . . ,x^) . Suppose the covariance matrix E is 
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( 4 ) 



The 



variance of f (Xj^| . . . ,x^) is given by: 

^k+1’ • • • »^n^ " ^11 “ ^12 ^22 ^21 



To simplify this expression, the following theorem is used 

[Ref. 4]: |e| = j^^i “ ^12 ^22^ ^21^ ’^^22^* square 

matrix ^^22 nonsingular since the covariance matrix E is 
positive definite. The variance (M can be written as: 



^(k) 

^^\l^k+l’ • • • ’^n^ " I j(k+l) 



(5) 



The mean of f (Xj^| Xj^_j_^, . . . ,x^) is given by 



^k+1’ • • • ’^n^ ^12 ^22 



Direct substitution of components of the partitioned covari- 
ance matrix (3) into the latter expression yields: 



^k+1' • • • ^^k,k+l ^k,k+2 •'* ^k,n^^^ 



(k+lX ^ nr -I 
^k+1 



X 



k+2 



X 

n 



A simplification of this expression can be made as follows. 
First, examine (E^^'*’^^) . By using standara matrix 

algebra, it can be shown that: 
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Next, perforin the matrix multiplication 
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a ] [Z^^'*'^^] 
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-1 



This yields a vector of dimension lx (n - k + 1). The 

transpose of this vector is the product of the factor 
1 
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Multiply the vector (6) by 
results are obtained: 



, and the following 



E(X, I X, ,x ) = 

k' k+1’ ’ n 






( k + 1 ) I J ^ ^1 ^ ^ k +1 ^ ^2 ^ \+2 



. . + (a^)x^] 



where (a^) indicates the first element of the vector (6), 
ia^) indicates the second element, etc. Upon close inspec- 
tion, it is apparent that (a^) is equivalent to k‘ 

In the same manner, (a^) is “^^+2 k’ **•» (a^) is 

_j;(k) mean of f(x, lx, x ) can thus be represented 

n,k k' k+1’ ’ n ^ 



as : 
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(7) 



The use of these expressions in generating the components 
of X is now described. Consider the, transformation 
Xj^ “ ^k ^ ^k’ ^ ~ 1>2,... n, where the Y's are indepen- 
dent normal variates with zero mean and unit variance. 

Using expressions (5) and (7) for and respectively, 
this transformation may be given explicitly as: 



X = Y /d 
n n n^n 



( 8 ) 
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(9) 



for k = l,...,n-l. Use of the equations (8) and (9) there- 
fore provides a method of generating samples of a random 
vector X, using univariate generators for the Yj^’s. 



C. TRIANGULARIZATION METHOD 



The covariance matrix, designated as E, is a symmetric 
positive definite matrix. It can be written as E = TT' 
where T is a lower triangular matrix [Ref. 6] , 
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X is obtained by making the non-singular transformation 

X = TY where Y Is distributed N(0,I). By previous remarks, 

X is distributed N(£,TT'). The matrix T can be obtained 

from E by using the so called "square root method" [Ref. 6]. 

Equations for t. . in terms of o. . are as follows: 
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The vector Y consists of independent normal components with 
zero mean and unit variance. These may be generated using 
a univariate normal random number generator. Standard 
matrix multiplication of T and samples of Y yields samples 
of the random vector X. 
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III. METHOD EVALUATION 



Computer runs were obtained using each of the three 
methods, with various size covariance matrices. Data con- 
cerning these runs is compiled in Tables I, II, and III. 

A brief description of the procedure employed to obtain 
execution time and memory space requirements for the tables 
follows. The computer programs can be subdivided into three 
basic portions. The first portion consists of the data 
that must be entered into the computer. The second portion 
consists of all required steps, prior to the use of the 
basic univariate random number generator, necessary to 
implement one of the three methods described in this study. 
The third portion consists of the standard Gaussian random 
number generator and those, steps necessary for generating 
samples of the random vector, X, such that a matrix is 
filled with one thousand random vectors. Two times were 
considered when evaluating the programs: the "set-up" time 
and the "repetition" time. These times correspond to 
execution time for the second and the third portion respec- 
tively of the computer program. Similarly, the memory 
space requirements shown in the tables make use of two 
numbers. The first number provides an indication of the 
amount of space required for the computer programs for 
each of the three methods. This number includes the sub- 
routine EIGEN in the rotation method and the function GRN 
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in all cases. The second number Is the total space 
required for the program plus any external function used 
such as square roots, absolute values, exponentials, input 
and output devices, etc. The second number varies as to 
the computer system being used and in this case the computer 
was the IBM 36 O/ 67 . The headings on two columns of Tables I, 
II, and III are labeled "Specific Conditioning" and 
"Specific Trlangularlzatlon" respectively. These refer to 
computer programs written specifically for 2x2 and 3x3 
covariance matrices, and are adaptations of the general 
programs . 

Based on the data contained in Tables I, II, and III, 
the best method to use for generating random vectors from 
covariance matrices dimensioned greater than 3x3 appears 
to be the trlangularlzatlon method. This method requires 
less memory space and considerably less execution time than 
the other two methods. If the user is interested in 
generating samples from a bivariate or trivarlate normal 
distribution, then either a specific conditioning or triang- 
ularization method can be used. Both methods require 
about the same amount of execution time and their space 
requirements are essentially the same. 

All methods use, as a basis, the univariate normal 
random number generator. Therefore, in order to establish 
a lower bound on times required in generating normal 
random vectors, the time required to generate normal random 
numbers was obtained. These numbers should serve as lower 
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bounds on the times required to generate a sample of one 
thousand repetitions of a 10, 5, 3> and 2 element vector 
respectively. The lower bounds obtained in this way are: 

10,000 numbers - 1.067^0 seconds 

5.000 numbers - 0.53209 seconds 

3.000 numbers - 0.32302 seconds 

2.000 numbers - 0.20332 seconds 

A specialized covariance matrix, (identity matrix), 
was used as an input to all three methods to provide the 
least cumbersome, mathematically speaking, of all matrices 
that would be encountered. The results of this test are 
tabulated in Table III. Each method maintained its overall 
ranking with respect to memory space and time requirements. 
The rotation method displayed a sharp decrease in set-up 
time which is reasonable as eigen values are easier (faster) 
to compute in this case. In general, each of the three 
methods required essentially the same amount of time for 
1000 repetitions using a typical covariance matrix as for 
the specialized covariance matrix. 
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TABLE I 



MEMORY SPACE REQUIREMENTS 

The first number is the amount of core space required 
for the computer program and the second number is the 
total amount of core space required in execution of the 
program. This second number includes, for example, 
external functions and input-output devices. Each 
space requirement is in bytes. 



Matrix 

Size 


Conditioning 


Trlangular- 

ization 


Rotation 


2 x 2 


11 , 792 / 37,168 


9 , 528 / 33 , 86^1 


11,952/36,288 


3 x 3 


15,880/^1,256 


13,568/37,90^ 


16,008/40,344 


5 x 5 


2it, 160/^9, 536 


21 , 712 / 46 , 0^8 


24 , 192 / 48,528 


10 X 10 


i^ 5 ,^ 2 i|/ 70,800 


42,336/66,672 


44,992/69,328 










- _ - 


Specific Cond. 


Specific Trlang. 


- - - - 


2 x 2 


8,752/33,088 


8,744/33,080 


----- 


3 x 3 


12,976/37,312 


12,960/37,296 


- - - - 
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TABLE II 



EXECUTION TIME (TYPICAL COVARIANCE MATRIX) 

The first number Is the set-up time, and the second 
number is the repetition time. Each execution time is 
in seconds. 



Matrix 

Size 


Conditioning 


Triangular- 

ization 


Rotation 


2x2 


. 00273 /. 33985 


.00106/. 33800 


.00496/. 37666 


3 x 3 


.007^1/. 55383 


. 00145/. 52411 


.01398/. 62023 


5 x 5 


.03975/1.06335 


. 00247 /. 99554 


.06117/1.23357 


10 X 10 


.70306/2.75560 


.00795/2.52280 


.41652/3.72289 










- 


Specific Cond. 


Specific Triang. 


----- 


2x2 


. 00091/. 252il6 


. 00101/. 27487 


- - - - 


3 x 3 


.00101/. 39530 


.00109/. 39572 


- - - - 
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TABLE III 



EXECUTION TIME (SPECIALIZED COVARIANCE MATRIX I) 

The first number is the set-up time, and the second 
number is the repetition time. Each execution time is 
in seconds . 



Matrix 

Size 


Conditioning 


Triangular- 

ization 


Rotation 


2x2 


.00278/. 3^299 


.00122 /. 33509 


. 00171/. 36376 


3x3 


.00793/. 55869 


. 00161 /. 5^223 


. 0023 V . 60250 


5x5 


. 0321 ^/ 1.06025 


. 00260 /. 98^1 il 9 


.00403/1.195^5 


10 X 10 


. 33795 / 2.69337 


.00733/2.^126^5 


. 01162 / 3.41616 










^ ^ 


Specific Cond. 


Specific Triang. 


^ - 


2x2 


.oolov . 25215 


. 00078 /. 2 i <921 


- - - - 


3x3 


. 00135/. 37903 


. 00096/. 38795 


- - - “ 
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APPENDIX A 



DISCUSSION OF THE COMPUTER PROGRAMS 

All three computer programs were written In Fortran IV. 
They are not optimal, i.e.. Improvements to the programs 
are possible. The results of these Improvements could 
change the value of memory space requirements and program 
execution times, but should not change the comparative 
ranking among the three methods. 

Two univariate normal random number generators that 
were available for use in the present study were GAUSS and 
GRN. The subroutine GAUSS is part of the IBM scientific 
subroutine package and is based on the mean value theorem. 
The function GRN was compiled by the Naval Postgraduate 
School staff, and is based^ on the G. Marsaglia technique 
[Ref. 7] . The execution time of GRN is approximately ten 
times faster than GAUSS, and for this reason it was used 
in all the programs. The rotation method uses the sub- 
routine EIGEN which is also available as part of the IBM 
scientific subroutine package. Three subroutines, (DTERM, 
REDMAT, and COFACT), were developed for use in the con- 
ditioning method. The subroutine DTERM computes the deter- 
minant of an n-square matrix. The subroutine COFACT 
removes the elements from the 1^^ row and column of an 

n-square matrix (Z), and provides the (n-l)-square matrix 
that is required in the computation of the cofactor Z. .. 
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The subroutine REDMAT reduces a matrix to the desired size. 

It Is capable of providing the user with a from an 

n X n dimension matrix where k £ n. 

Each method was evaluated using four different size 
matrices as listed in Tables I, II, and III. These were 
selected as being typical of the range of matrix sizes that 
might be encountered In practice. However, the programs 
were written In general terms so there are no program 
imposed restrictions on matrix size. 

As mentioned previously, specific computer programs 
were written for the 2x2 and the 3x3 covariance matrices 
using the conditioning and trlangularizatlon method. This 
was done in an attempt to evaluate the possible reduction 
In memory space and execution time under those obtained 
using the general programs with Input dimensions of 2 and 
3. Of the three methods considered, these two appeared to 
be adaptable to simpler, more concise programs for these 
small dimensions. Writing specific programs for covariance 
matrices dimensioned greater than 3x3 would become dif- 
ficult, and a general program would probably have to be 
used. No attempt was made to write a specific program for 
the rotation method since computation of eigen values and 
vectors are involved, and there were no apparent means of 
reducing the computation time under that of using the general 
subroutine EIGEN. 
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The computer programs, used in generating multivariate 



normal random 
Included as a 



vectors for each of the three methods, 
portion of this study. 



2^ 



are 



COMPUTER PROGRAMS 



THIS PROGRAM UTILTZFS THE ROTATION METHOD TO GFNERATP 
MULTIVARIATE NORMAL RANDOM VECTORS. 



DIMENSION 7(3,3), P(3,3), XA(3,1000), R(5), Y(3), T(3) 
READ (5,1) N 

1 FORMAT (15) 

m 2 I = 1,N 

2 READ (5,3) (7( I , J) , J = 1 ,N) 

3 FORMAT (3F8.0) 

MV = 0 

K = 1 

DO 5 J = 1,N 
DO 4 I = 1,J 
B(K) = Z( I ,J) 

K = K + 1 

5 CONTINUE 
LB = K - 1 

CALL EIGEN (B,R,N,MV) 

MA = 1 
JA = 1 

DO 5 LA = 1,N 
Y( JA ) = R( MA ) 

Y( JA) = SQRT(Y( JA) ) 

JA = JA + 1 

6 MA = HA + JA 

IDUMMY = 0 

DO 10 MT = 1,1000 

DO 7 JP = 1,N 

X = ( GRN( IDUMMY) ) *Y( JP) 

7 T( JP ) = X 

DO 9 I = 1,N 
SUM = C.O 
DO 8 J = 1,N 

8 SUM = SUM R( I , J) Y T( J) 

9 XA( I,MT ) = SUM 

10 CONTINUE 
STOP 
END 
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THIS PROGRAM UTILIZES THE TP I ANG' IL A P T Z AT I PM MPTHOP 
TO GENERATE MULTIVARIATE NORMAL RANDOM VECTORS. 



niMENSICN Z(3,3), C(3,3), XV(1000t3). T(3) 
READ ( E,1 ) N 

1 FORMAT (18) 

READ (5,2) ((7(I,J), J = 1,N), I = 1, N) 

2 FORMAT (3F8.0) 

DO A. I = l,N 

C( I , 1 ) = 7(1,1) /SORT (7(1,1)) 

IF ( I +1 .GT. N) GO TO 
JE = I + 1 
no 3 JR = JP,N 
C( I, JR ) = 0.0 

3 CQNTINLF 

A CONTINUE 

DO 10 I = 2,N 
IF ( T - 3) 8,5,5 
5 JT = I - 1 

DO 7 J = 2,JT 

SUM = 0.0 
JD = J - 1 
DO 6 KA = i,jn 

SUN = SUN + C(I,KA) * C(J,KA) 

5 CONTINUE 

C( I , J ) = ( Z( I , J) - SUN) /c ( J , J) 

7 CONTINUE 

8 SUM = 0.0 

JF = I - 1 

DO 9 K = 1,JF 

SUM = SLM + C(I,K) C(I,K) 

9 CONTINUE 

C( I ,I ) = SORT (7(1,1) - SUM) 

10 CONTINUE 
IDUMMY = 0 

DO lA JP = 1,1000 

DO 11 JC = 1 ,N 

X = GRN(IDUMMY) 

11 T( JC ) = X 

DO 13 IT = 1,N 
SUMM = 0.0 
DO 12 J = 1,IT 

12 SUMM = SUMM + C(IT,J) ^i-TIJ) 

XV(JP,IT) = SUMM 

13 CONTINUE 

14 CONTINUE 
STOP 

END 
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THIS PROGPiM UTILIZES TH<= CDNOITIOHAL MFTHOD TO 
GENERATE MULTIVARIATE NORMAL RANDOM VECTORS. 



DIMENSION Z(3,3), A(?,3), R(3,3), Y(3), R(3), T(3)t 
1XA(1000,3), ZJ(3,3) 

READ (5,1) N 

1 FORMAT (lA) 

READ (5,2) ((Z(I,J), J = 1,N), I = 1,N) 

2 FORMAT (3F8.0) 

M = N - 1 

DO 6 L = 1,M 

K = N - L 

LT = N - K + I 

CALL REDMAT (Z,N,LT,K,A) 

IR = K + 1 
DO 3 IV = TR, N 

MV = LT - 1 

IF = IV - K + 1 

IG = 1 

CALL COFACT ( L T , A , B , MV , I F , I G ) 

CALL DTEPM (NV, B, D, NV) 

IT = IF + IG 

3 ZJ(K,IV) = ((-l)’!'*IT) Jf D 

CALL OTFPM (LT,A,D,LT) 

DNUM = D 

DO A IS = 1,N 
DO A JS = 1,N 
A( IS, JS) = 0.0 

B( IS, JS ) = 0.0 

A CONTINUE 

CALL REDMAT ( Z ,('! ,NV, IR ,A) 

CALL DTEPM(NV,A,D,NV) 

T(K) = D 
DO 5 IC = 1,N 
DO 5 JC = 1,N 
A( IC , JC ) = 0.0 

5 CONTINUE 

6 R(K) = SORT(DNUM/T( K) ) 

AO = SORT! Z( N,N) ) 
lOUMMY = 0 

DO 10 NOP = 1,1000 
DO 7 NO = 1,N 
X = GRN(IDUMMY) 

7 Y(NO) = X 
XA(NOP,N) = AD * YIN) 

LA = M - 1 

DO 9 NORM = 1,LA 
SUMN = 0.0 
K = N - NORM 
LEE = K + 1 
DO 8 NSYL = LEE ,N 

8 SUMN = ZJ(K,NSYL) XA( NOP, NSYL) + StJMN 

9 XA(NOP,K) = R(K) A Y(K) - (SUMNZT(K)) 

10 CONTINUE 
STOP 

END 
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SUBPOUTINF DTFPM (N,A,D,M) 
niMPMSlOM A(M»M) 

DO = 1.0 
on lo L = ItN 
KP = 0 
7 = 0.0 
DO 12 K = L,N 

IF (Z - ARS( A(K,L) ) ) 11,12,12 
Z = ABS(A(K,D) 

KP = K 
CONTINUF 

IF (L - KP) 13,15,15 
DO lA J = L,M 
Z = A( L , J ) 

A(L,JI = A(KP,J) 

A( KP , J ) = Z 
DO = - DO 

IF (L - N) 16,20,20 
LPl = L + 1 
DO 19 K = LPl , M 
IF (A(K,L) ) 17,19,17 
RATIO = A(K,L) /A( L,L) 
no 18 J = LPl ,M 

A(K,J) = A(K,J) - RATIO >!' A(L,J) 
CONTINUF 
DO 21 K = 1,N 
DD = DD « A( K ,K ) 

D = DD 
RETURN 
END 



SUBROUTINE REDMAT ( 7 , N , L 7 , 1 B , A ) 
DIMENSION Z(N,N), A(L7,LZ) 

I = 0 

DO 23 L = 1,N 
IF (L .LT. IB) GO TO 23 
1 = 1 + 1 
J = 0 

DO 22 K = 1,N 

IF ( K .LT. IB) GO TO 22 

J = J + 1 

A( I , J ) = Z( L ,K ) 

CONTINUF 

CONTINUF 

RETURN 

END 



SUBROUTINE COFACT ( NA , Z , B , NV , I K , I L ) 
DIMENSION Z(NA,NA), B(NV,NV) 

I = 0 

DO 25 K = 1,NA 
IF (K .EO. IK) GO TO 25 
1 = 1 + 1 
J = 0 

DO 2*- L = 1,NA 

IF (L .EO. ID GO TO 24 

J = J + 1 

B( I ,J) = 7( K,L) 

CONTINUE 

CONTINUE 

RETURN 

END 






28 



BIBLIOGRAPHY 



1. Wold, H., Tracts for Computers , No. XXV, Cambridge 

Press, 1955. 

2. Scheuer, E. M. and Stoller, D. S., "On the Generation 

of Normal Random Vectors," Technometrics , v. ■4, 
no. 2, p. 278-281, May 1962. 

3. Muller, M. E., "A Comparison of Methods for Generating 

Normal Deviates on Digital Computers," Journal of 
the Association for Computing Machinery"] V. 6 , no. 3j 
p. 376-383, July 1959. 

Anderson, T. W., An Introduction to Multivariate 

Statistical Analysis , p. 17-19> App. A, Wiley, i960. 

5. Graybill, F. A., An Introduction to Linear Statistical 

Models , p . 5 , v"! I, McGraw-Hill, 196I . 

6. Graybill, P. A., Introduction to Matrices with Applica - 

tions in Statistics , p. 298, Wadsworth, 1969 • 

7. Marsaglia, G., "A Past Procedure for Generating Normal 

Random Variables," Communications of the ACM, v. 7> 
no. 1, p. 4-10, J anuary 1964 . 

8. Howe, J. E. , The Generation of Random Numbers from 

Various Probability Distributions , Masters Thesis, 
Naval Postgraduate School, Monterey, I965. 



29 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Documentation Center 2 

Cameron Station 

Alexandria, Virginia 2231^ 

2. Library, Code 0212 2 

Naval Postgraduate School 

Monterey, California 939^0 

3. Department of Operations Analysis, Code 55 1 

Naval Postgraduate School 

Monterey, California 939^0 

Associate Professor D. R. Barr, Code 55 bn 1 

Department of Operations Analysis 
Naval Postgraduate School 
Monterey, California 939^0 

5. LCDR N. L. Slezak, USN 1 

Box 3^ 

Milligan, Nebraska 68^06 



30 



UNCLASSIFIED 

Socuhlv Classification 



DOCUMENT CONTROL DATA - R & D 

(Security classitication of ttttc, body of abstract and indexinfi annotation rtfust be entered when the overall report Is c lassified) 



i. originating activity (Corporate author) 


2a. REPORT security CLASSIFICATION 


Naval Postgraduate School 


Unclassified 


2b, GROUP 


Monterey, California 939^0 





3 REPORT TITLE 



A COMPARISON OP METHODS FOR GENERATING 
MULTIVARIATE NORMAL RANDOM VECTORS 



4 DESCRIPTIVE NOTES (Type of report and, inclusi ve dates) 

Master’s Thesis: September 1Q70 

5 . AU T HO R IS > fF/rsf name, middle initial, last name) 

Norman Lee Slezak, Lieutenant Commander, United States Navy 



6 REPORT DATE 


7a. TOTAL NO. OF PAGES 


7b. NO, OF REFS 


September 1970 


30 


00 


ea. CONTRACT OR GRANT NO- 


9a, ORIGINATOR'S REPORT NUMBER(S) 


b. P ROJEC T NO. 






C. 


9b. OTHER REPORT NO(S) (Any other nuonbere that may be assigned 
this reporf^ 


d. 







10. DISTRIBUTION STATEMENT 



This document has been approved for public release and sale; 
Its distribution Is unlimited. 



n. SUPPLEMENTARY NOTES 



12. SPONSORING MILITARY ACTIVITY 



Naval Postgraduate School 
Monterey, California 939^0 



Three methods for generating outcomes on multivariate 
normal random vectors are presented. A comparison Is made 
to determine which method requires the least computer 
execution time and memory space when utilizing the 
IBM 360 / 67 . All methods use as a basis a standard Gaussian 
random number generator. Results of the comparison study 
Indicate that the method based on triangular factorization 
of the covariance matrix generally requires less memory 
space and computer time than the other two methods . 




31 



UNCLASSIFIED 



Security Classification 



A-31408 



UNCLASSIFIED 

Security Classification 



KEY WORDS 



ROLE W T 



MULTIVARIATE NORMAL DISTRIBUTION 
RANDOM NUMBER GENERATOR 
NORMAL DISTRIBUTION 
MULTIVARIATE NORMAL GENERATOR 



3 .?olM473 (^ack) 



01 Oi *807-6821 



32 



TIMCLASSIFIED 



Security Classification 



A-31 409 



120094 




Thesi s 

S5709 

C.l 



7LH 



Slezak _ , ^ 

A comparison or 
„,ethods for general. nc 
multivariate normal ^ 
random vecto^2 5 15 

V r 






Thesis 120094 

S5709 Slezak 

C.l A compari son of 

methods for generating 
multivariate normal 
random vectors. 



itu'sSS '0^» 

A comp.Hision ot methoiis tor generjting 




3 2768 001 00617 4 

DUDLEY KNOX LIBRARY 



